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ABSTRACT 

According to the no-hair theorem, astrophysical black holes are uniquely described by their masses 
and spins. An observational test of the no-hair theorem can be performed by measuring at least three 
different multipole moments of the spacetime of a black hole and verifying whether their values are 
consistent with the unique combinations of the Kerr solution. In this paper, we study quasi-periodic 
variability observed in the emission from black holes across the electromagnetic spectrum as a test of 
the no-hair theorem. We derive expressions for the Keplerian and epicyclic frequencies in a quasi-Kerr 
spacetime, in which the quadrupole moment is a free parameter in addition to mass and spin. We show 
that, for moderate spins, the Keplerian frequency is practically independent of small deviations of the 
quadrupole moment from the Kerr value, while the epicyclic frequencies exhibit significant variations. 
We apply this framework to quasi-periodic oscillations in black-hole X-ray binaries in two different 
scenarios. In the case that a pair of quasi-periodic oscillations can be identified as the fundamental 
g- and c-modes in the accretion disk, we show that the no-hair theorem can be tested in conjunction 
with an independent mass measurement. If pairs of oscillations are identified with non-parametric 
resonance of dynamical frequencies in the accretion disk, then testing the no-hair theorem also requires 
an independent measurement of the black-hole spin. In addition, we argue that VLBI observations of 
Sgr A* may test the no-hair theorem through a combination of imaging observations and the detection 
of quasi-periodic variability. 

Subject headings: accretion, accretion disks — black hole physics — gravitation — hydrodynamics — 
stars: individual (Sgr A*) — X-rays: binaries 
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1. INTRODUCTION 

In general relativity, the no-hair theorem rests on the 
assumptions that any spacetime singularity must be en- 
closed by an event horizon (the cosmic censorship con- 
jecture, Penrose 1969) and that the exterior spacetime is 
free of closed timelike curves. 

If these assumptions are satisfied, astrophysical black 
holes are uniquely characterized by their mass M and 
spin J and are described by the Kerr metric (Israel 1967, 
1968; Carter 1971, 1973; Hawking 1972; Robinson 1975). 
Mass and spin are the first two multipole moments of a 
black- hole spacetime, and all higher order moments are 
fully specified by the relation (Hansen 1974) 



Mi+iSi = M(ia) 1 . 



(1) 



Here, a = J/M is the spin parameter, and the multipole 
moments are written as a set of mass multipole moments 
Mi, which are nonzero for even values of I, and a set 
of current multipole moments Si, which are nonzero for 
odd values of I. Since the multipole moments of a Kerr 
black hole have to be related by expression (fTJ), the no- 
hair theorem can be tested by measuring (at least) three 
such moments (Ryan 1995). 

Potential tests of the no-hair theorem have been sug- 
gested using observations of gravitational waves from ex- 
treme mass-ratio inspirals (Ryan 1995, 1997a, 1997b; 
Barack & Cutler 2004, 2007; Collins & Hughes 2004; 
Glampedakis & Babak 2006; Gair et al. 2008; Li & 
Lovelace 2008; Apostolatos et al. 2009; Vigeland & 
Hughes 2010), from electromagnetic observations of stars 



in close orbits around Sgr A* (Will 2008; Merritt et al. 
2010), as well as of pulsar black-hole binaries (Wex & 
Kopeikin 1999). 

In the first part of this series of papers (Johannsen 
& Psaltis 2010a, hereafter Paper I), we investigated a 
framework for measuring three multipole moments of 
black holes with observations in the electromagnetic 
spectrum. We used a quasi-Kerr metric (Glampedakis & 
Babak 2006), which contains an independent quadrupole 
moment of the form 



Q = -M i 



eM 2 ) 



(2) 



where a potential deviation from the Kerr metric 
quadrupole is parameterized in terms of e. This param- 
eterization includes the Kerr metric as the special case 
e = 0. 

If a measurement yields e ^ 0, then the compact ob- 
ject cannot be a Kerr black hole. Within the framework 
of general relativity, it can only be a different type of 
star, a naked singularity, or else an exotic configuration 
of matter (see Collins & Hughes 2004; Hughes 2006). If, 
however, the compact object is otherwise known to pos- 
sess an event horizon and a regular spacetime, then a 
quadrupole deviation implies that the no-hair theorem is 
incorrect and that general relativity itself breaks down 
in the strong-field regime. We refer to black holes of 
this kind as "quasi-Kerr black holes" (Paper I; see, also, 
"bumpy black holes" Collins & Hughes 2004; Vigeland 
k Hughes 2010). 

As a first astrophysical application of our framework 
(Johannsen & Psaltis 2010b, Paper II), we simulated im- 
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ages of such objects and demonstrated their dependence 
on the quadrupolc moments of their spacetimes. In par- 
ticular, we showed that a bright and narrow ring outside 
of the shadow of the black hole (hereafter the ring of light; 
see Johannsen & Psaltis 2010c) with a diameter of about 
10M (see also, e.g., Beckwith & Done 2005) has a shape 
that depends uniquely on the mass, spin, quadrupole mo- 
ment, and inclination of the black hole. In this paper, 
we show how quasi-periodic variability observed across 
the electromagnetic spectrum in the emission from black 
holes can be used to test the no-hair theorem. 

Quasi-periodic oscillations (QPOs) have been observed 
in several galactic binaries with the Rossi X-ray Timing 
Explorer (RXTE; see Remillard & McClintock 2006; see 
van der Klis 2006 for definitions and analysis techniques) 
and in Active Galactic Nuclei (AGN) with XMM- Newton 
(Gierlihski et al. 2008). While potential QPO signals 
in AGN are often obscured by red noise (Benlloch et 
al. 2001; Vaughan & Uttley 2005, 2006), galactic black 
holes usually reveal much cleaner signals. In the case 
of galactic black holes, QPOs arc transient phenomena 
that occur during mostly nonthermal states of the black- 
hole accretion disk and during state transitions; they fall 
into two general classes: high-frequency (about 40-450 
Hz) QPOs and low-frequency (about 0.1-30 Hz) QPOs 
(Remillard & McClintock 2006). 

The physical origin of the observed QPOs is not well 
understood. In a thoroughly developed hydrodynamic 
model, they attributed to normal modes of oscillation 
trapped by general-relativistic effects in the accretion 
disks around black holes (see Wagoner 1999 and Kato 
2001 for reviews and references therein). Expressions of 
these modes and the corresponding oscillation frequen- 
cies have been derived for the case of modified Newto- 
nian potentials (Kato & Fukue 1980; Okazaki et al. 1987; 
Kato 1990; Nowak & Wagoner 1991, 1992, 1993) and in 
full general relativity (Perez et al. 1997; Silbergleit et 
al. 2001; Ortega-Rodriguez et al. 2002). Since these fre- 
quencies depend primarily on the mass and the spin of 
the black hole and only very little on the speed of sound, 
quasi-periodic oscillations provide a laboratory for tests 
of general relativity (see Psaltis 2003, 2008 for reviews). 

Of special interest are gravity modes in the equato- 
rial plane, the so-called (/-modes (Perez et al. 1997; see, 
however, Li et al. 2003), which are trapped near the 
inner edge of the accretion disk, as well as corrugation 
modes or c-modes (e.g., Silbergleit et al. 2001), that pre- 
cess slowly around the angular momentum vector of the 
black hole. Both of these modes lead to a modulation 
of black-hole X-ray spectra; g-modes usually cover the 
largest area of the disk near the temperature maximum, 
while c-modes affect the projected area of the disk (see, 
e.g., Wagoner 1999). Both the g- and the c-modes are 
related directly to the epicyclic frequencies of particles 
on (nearly) circular equatorial orbits (Perez et al. 1997; 
Silbergleit et al. 2001). The most robust and observable 
mode is expected to be the axisymmetric g-modc, which 
has also been seen in various hydrodynamic simulations 
(Reynolds & Miller 2009; Mao et al. 2009; Chan 2009; 
see, also, Wagoner 2008). 

In a different approach, since QPOs have been observed 
in pairs with frequency ratios of w 3/2 in several sources 
(see Remillard & McClintock 2006), they have also been 
modeled as nonlinear resonances among the Keplerian 



and epicyclic frequencies (Kluzniak & Abramowicz 2001; 
Abramowicz et al. 2003). The frequencies predicted by 
both of these models are consistent with observations and 
have been used in each model to constrain the spin of 
the black-hole X-ray binary GRO 1655-40 (Abramowicz 
& Kluzniak 2001; Wagoner et al. 2001). 

Variability has also been observed in the emission of 
Sgr A* in the radio, millimeter, NIR, and X-ray bands 
with timescales ranging from minutes to hours (e.g., 
Baganoff et al. 2001; Aschenbach et al. 2004; Genzel et 
al. 2003; Ghez et al. 2004; Belanger et al. 2006; Meyer et 
al. 2006; Yusef-Zadeh et al. 2006; Marrone et al. 2006; 
Hornstcin et al. 2007). Possible explanations of such 
variability include models of density inhomogencitics 
("hot spots") orbiting around Sgr A*. Doeleman et 
al. (2009) investigated the prospects of (sub— ) millimeter 
very-long baseline interferometry (VLBI) to detect peri- 
odicity in the emission of Sgr A* and to measure its spin. 

In this paper, we derive expressions for the Keplerian 
and epicyclic frequencies of circular equatorial motion in 
the quasi-Kerr spacetime. We discuss the properties of 
these frequencies and analyze their distinct dependencies 
on the mass, spin, and quadrupole moment of the black- 
hole spacetime. In particular, we demonstrate how this 
formalism can be applied in order to test the no-hair the- 
orem with the quasi-periodic variability observed from 
galactic black holes, AGN, and Sgr A*. 

In Section 2, we derive expressions for the Keplerian 
frequency and the radial and vertical epicyclic frequen- 
cies in quasi-Kerr spacetimes and analyze the stability of 
circular orbits. In Section 3, we apply our expressions to 
quasi-periodic oscillations in black-hole X-ray binaries. 
We formulate our conclusions in Section 4. 

2. DYNAMICAL FREQUENCIES IN QUASI-KERR 
SPACETIME 

In this section, we systematically derive expressions 
for the Keplerian frequency fi^ of a particle on a circular 
equatorial orbit around a quasi-Kerr black hole as well 
as for the frequencies of small oscillations k t and fig in 
the two directions perpendicular to such an orbit. For 
an alternative computation of the frequencies 17 ^ and K r , 
see Glampcdakis & Babak (2006). 

Our starting point is the Kerr metric g^ h , which 
in Boyer-Lindquist coordinates takes the form (e.g., 
Bardeen, Press, & Teukolsky 1972) 



ds 2 



1 



2Mr 



, 2 f 4Mar sin 2 0\ 
v / '" ( v ) (Ud<> 



with 



— dr 2 +Y,d0 2 + [r 2 + a 2 + 



2M a 2 r sin 2 6 



sin 



A = r 2 - 2Mr + a 2 , 
£ = r 2 + a 2 cos 2 6. 



(3) 



(4) 



The quasi-Kerr metric g^ h in Boyer-Lindquist coordi- 
nates is given by (Glampcdakis & Babak 2006) 

ffa Q b K =.9a K b+^ab. (5) 

In contravariant form, h ah is 

h u = (1 - 2M/r)- 1 [(1 -3 cos 2 (9) Ji(r)] , 



TESTING THE NO-HAIR THEOREM: III. QUASI-PERIODIC VARIABILITY 

h" = (1 - 2M/r) [(1-3 cos 2 6) T x (r)] , 



[(l-3cos 2 0) J- 2 (r)] , 



9 ■ 2 

r z sin 



-[(l-3cos 2 0)j- 2 (r)], 

= 0. (6) 



The functions Fi^fr) are given in Appendix A of 
Glampcdakis & Babak (2006). In its full form, the quasi- 
Kerr metric is valid only for slowly rotating black holes 
with values of the spin a < 0AM. Note, however, that 
the unperturbed spacetime (i.e., e = 0) is formally cor- 
rect up to the maximum value of the spin. 

Since the quasi-Kerr metric is stationary and axisym- 
metric, particle trajectories in this metric are character- 
ized by three integrals of motion. For a particle with 
4-momentum 

dx a 



P 



/'■ 



dr ' 



(7) 



these are its rest mass /i (which we will set equal to unity 
from here on), energy E — —pt, and angular momentum 
about the 2- axis L z = p^. 

We use the conservation of energy and axial angular 
momentum to express the momentum components p t and 
p^ of a particle in the form 



P 



P f 



g<p<pE + g t< f>L z 
9ttg<p<p - g 2 ^ 

g t< pE + g tt L z 



(8) 



(9) 



We then bring the conservation equation for the parti- 
cle's rest mass, 

-1, (10) 



9^p v 



into the form 



2 ( 5rr 



= ^ [-9tt(pT ~ 2 W V - 9M(p*Y -1\= V cS (11) 

with Vcff playing the role of an effective potential for the 
particle motion in the coordinates r and 9. 

The radial motion in the equatorial plane can be ana- 
lyzed in terms of the equation 



1 (dr 

2 \ dt 



t\2 - V eSi 



g™(p 



(12) 



while the 0-motion (for constant radii r) is governed by 
the equation 



1 (dQ_ 

2 [dt 



Vol 



,t^2 eff- 



geeip 

In Paper I, we solved the equations 



and 



off 



dr 



(13) 

(14) 
(15) 



which characterize circular equatorial orbits, and ob- 
tained the energy 



E 



r 3 ' 2 - 2Mr 1 / 2 ± aM 1 / 2 
r 3/4^ r 3/2 _ 3Mr i/2 ± 2aM 1 / 2 

he 



32M 2 r 3 /2( r _3M) 3 /2 
2M(6M 4 + 14M 3 r - 41M 2 r 2 + 27Mr 3 - 6r 4 ) 

r 2 (6r 3 - 33Mr 2 + 66M 2 r - 48M 3 ) In 



r-2M 



and axial angular momentum 

M 1 /2( r 2 T 2aM 1 /V/ 2 + a 2 ) 



(16) 



U = ±- 



" r 3 / 4 Vr 3 / 2 - 3MrV2 ± 2a M 1 l 2 
5e 



T 32M 5 / 2 (r-3M) 3 / 2 
2M(6M 4 - 7Af 3 r - 16M 2 r 2 + 12A/r 3 - 3r 4 ) 

-3r(r 4 - 5Mr 3 + 9A/V - 2Af 3 r - 6M 4 ) In 



r-2M 

These expressions are expansions to linear order in the 
parameter e, where we have neglected terms of the order 
ea. An expansion of this type is implicitly understood for 
all expressions throughout the paper in accordance with 
the form of the metric specified by expression ([5]) . Note, 
however, that we do not expand in the spin parameter a 
so that our results are correct for arbitrary values of the 
spin in the special case e = 0. 

2.1 THE KEPLERIAN FREQUENCY 

For a particle moving on a circular orbit in the equato- 
rial plane, the time and azimuthal components of the 4- 
momentum take the form (Glampedakis & Babak 2006) 



"4 



E(r 



2Ma 



(aE - L x ) 



A 



2M 



(aE - L z ) + U 



(r)E, 



(18) 



(19) 



where 



/a(r) 



5(r - M) 
~8Mr(r-2M) 



(2M 2 + 6Mr -3r 2 ) 



15r(r - 2M) 



h 3 (r) 



16M 2 
5 

" 8Mr 



In 



r - 2M J ' 
(2M 2 - 3Mr- 3r 2 ) 



15 
16M 2 



^r 2 - 2M 2 ) In 



r-2M 



(20) 
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From these expressions, we obtain the Keplerian fre- 
quency 

o _ # 

at p t 



±- 



5e 



I 2 ± aVM 32M 5 / 2 (r - 2M)r 5 / 2 
2M(-6M 4 + 8M 3 r + 2M 2 r 2 + 3Mr 3 - 3r 4 ) 



+ 3r(r - 2A/)(r 3 - 2M 3 ) In 



r-2M 



(21) 




0.02 



radius (M) 



Fig. 1. — Keplerian frequency f2^ of a particle on a circular equa- 
torial orbit around a quasi-Kerr black hole with a spin of a = 0.2M 
and for several values of the parameter e. The Keplerian frequency 
increases for decreasing values of the radius and for increasing val- 
ues of the parameter e. The dots mark the location of the ISCO, 
and the dashed lines correspond to unstable orbits. 

In Figure [T] we plot the Keplerian frequency as a func- 
tion of radius for several values of the parameter e and 
a spin of a — 0.2M . The frequency increases at smaller 
radii. The dependence on the quadrupolar parameter 
e is weak and only manifests at radii comparable to a 
few M. Increasing the value of the parameter e leads to 
larger Keplerian frequencies. 

2.2 THE RADIAL AND VERTICAL EPICYCLIC 
FREQUENCIES 

In order to calculate the radial and vertical epicyclic 
frequencies, we introduce small perturbations Sr and 59 
around a circular equatorial orbit located at radius rrj. 
We take the (coordinate) time derivative of equations 
([121) and (O) and obtain 



d 2 (Sr) 
dt 2 

d 2 (se) 
dt 2 



d 2 K 



cir 



dr 2 
d 2 V 



Sr. 



oil 



de 2 



so 



using the effective potentials V* s and VjL defined in equa- 
tions (jl"2j) and (fl~3|) . From these expressions, we derive 
the frequencies of small oscillations in the r and 9 direc- 
tions around a circular equatorial orbit as 



Or, 



dr 2 

de 2 



(24) 



(25) 



Here the second derivatives of the effective potentials are 
evaluated at the energy and axial angular momentum 
given by expressions (fl~6f and (flTj) . and for 9 = n/2. 
The result is: 



M (r 2 - 6Mr ± 8aVM - 3a : 
r 2( r 3/2 ±ay/M) 2 

5e 



16M 2 r 5 (r-2M) 



2M 



(48M 5 + 30M 4 r + 26M 3 r 2 - 127M V + 75Afr 4 - 12r 



+3r 2 (r - 2M) 2 (4r 2 - 13Mr - 2M 2 ) In 



2M 



and 



(26) 



M(r =F ^a\jMsJr + 3a 



r 2 (r 3 / 2 ±aVM) 2 
5c 



lQM 2 r 4 (r-2M) 
2M(6M 4 + 34M 3 r - 59M 2 r 2 + 33Mr 3 

+3r(2r-M)(r-2M) 3 ln 



6r 4 ) 



2M 



(27) 



To convert these expressions to cgs units, one needs to 
multiply by (c 3 /G) 2 , where c and G are the speed of 
light and the gravitational constant, respectively. In the 
special case where e = 0, these expressions coincide with 
the ones in Okazaki et al. (1987) and Kato (1990). Sim- 
ilar calculations have also been performed by Shibata 
& Sasaki (1998), Berti & Stergioulas (2004), Gair et al. 
(2008), and Sanabria-Gomez et al. (2010). 

We plot the radial epicyclic frequency n r as a function 
of radius for several values of the parameter e and a spin 
of a = 0.2M in FigureEJ The radial oscillation frequency 
increases with decreasing values of the radius and reaches 
a maximum at r ~ 7M for this spin. Decreasing values 
of the parameter e shift the maximum to smaller radii. 
In Figure [3] we plot the Lense-Thirring frequency 

n LT =fl^-Q e (28) 

as a function of radius for several values of the parame- 
ter e at a spin value a — 0.2M. The Lense-Thirring fre- 
quency describes the precession of the orbital plane of a 
particle moving around the black hole and vanishes in the 

(22) case °f a Schwarzschild black hole, where = The 
Lense-Thirring frequency increases with decreasing val- 
ues of the radius and of the parameter e. This frequency 

(23) depends significantly on the parameter e and changes 
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Fig. 2. — Radial cpicyclic frequency K T versus radius for a particle 
orbiting a quasi-Kcrr black hole with a spin of a = 0.2M and for 
several values of the parameter e. The radial epicyclic frequency 
increases with decreasing values of the radius and reaches a maxi- 
mum at r Kb 7M for this spin. Decreasing values of the parameter 
e shift the maximum to smaller radii. 



Fig. 4. — Location of the ISCO as a function of the spin a and 
the parameter e. The radius of the ISCO decreases with increasing 
values of the spin and decreasing values of the parameter e. The 
dashed line corresponds to the special case of the Kerr metric, while 
the shaded region is excluded and marks the part of the parameter 
space where higher order terms in the radial epicyclic frequency 
become important. 
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Fig. 3. — Lense-Thirring frequency f^LT = ^4> — &6 versus radius 
for a particle on an orbit around a quasi-Kerr black hole with a 
spin of a = 0.2M and for several values of the parameter e. The 
Lense-Thirring frequency increases with decreasing values of the 
radius and of the parameter e. The dots mark the location of the 
ISCO, and dashed lines correspond to unstable orbits. 

by a factor of ~ 5 at the ISCO for values of the pa- 
rameter e = ±0.2. Note that both the Lense-Thirring 
frequency and the radial oscillation frequency are signif- 
icantly smaller than the Keplerian frequency f2^. 



Fig. 5. — Radial epicyclic frequency function of radius 

for a black hole of spin a = 0AM and for several values of the 
parameter e. Zeros of the radial cpicyclic frequency correspond to 
marginally stable orbits. Small negative values of the parameter e 
lead to the emergence of a second marginally stable orbit in addi- 
tion to the ISCO. For values of the parameter e < —0.23, these two 
merge and all orbits are now stable against radial perturbations. 
The ISCO is then determined by the onset of a vertical instability. 

2.3 STABLE ORBITS IN QUASI-KERR SPACETIME 

In this section, we determine the location of the ISCO 
in the quasi-Kerr spacetime which separates the region 
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-0.4 -0.2 0.2 0.4 

spin (M) 



Fig. 6. — Maximal Keplerian frequency in Hz scaled appropri- 
ately for a lOM0-mass black hole as a function of the spin a and 
the parameter e. The maximal Keplerian frequency increases with 
increasing values of the spin and decreasing values of the parame- 
ter e. The dashed line corresponds to the special case of the Kerr 
metric. The shaded region marks the excluded part of the parame- 
ter space where higher order terms in the radial epicyclic frequency 
become important. 

of space where circular orbits are stable from the region 
where circular orbits are unstable. The location of this 
orbit, in turn, determines the maximum Keplerian fre- 
quency for that particular central object. 

In Paper I we found the location of the ISCO from 
equation (fT6| by solving 



In this paper, we determine the location of the ISCO by 
computing the radius r at which the frequency of small 
radial oscillations k t becomes imaginary, rendering cir- 
cular equatorial motion unstable (see, also, Shibata & 
Sasaki 1998; Berti & Stergioulas 2004; Gair et al. 2008). 

In Figured] we plot the location of the ISCO as a func- 
tion of the spin and the parameter e. The location of the 
ISCO shifts to smaller radii for increasing values of the 
spin as in the case of Kerr black holes, while the radius 
of the ISCO increases for increasing values of the param- 
eter e. We denoted the location of the ISCO in Figures[T] 
and [3] by dots. 

In addition to the region outside of the ISCO, we find 
a second branch of stable circular orbits similar to the 
ones reported in Gair et al. (2008). If the parameter 
e is sufficiently negative for a given value of the spin, 
we observe the emergence of this branch at radii smaller 
than the radius of the ISCO as well as of a second unsta- 
ble branch for which the epicyclic frequency fig becomes 
imaginary. 

We illustrate this behavior in Figure [5J where we plot 
the radial epicyclic frequency as a function of radius for 
a black hole with spin a = 0.4M and for several values 
of the parameter e. For a Kerr black hole, the radial 



epicyclic frequency has only one zero at a radius nsco ~ 
4.6M. If e < 0, the location of the ISCO decreases, while 
a second branch of stable circular orbits emerges in the 
neighborhood of the circular photon orbit. The radius 
of the second marginally stable circular orbit increases 
with decreasing values of the parameter e. For values of 
the parameter e < —0.23, the radial epicyclic frequency 
becomes strictly positive and these two marginally stable 
circular orbits disappear. In this case, all circular orbits 
are stable against radial perturbations, and the ISCO is 
determined by instability in the vertical direction. 

Gair et al. (2008) referred to the larger one of these 
additional two marginally stable circular orbits as OSCO, 
the outermost stable circular orbit. In all cases, however, 
r OSCO < r isco- Since a particle on a circular orbit in this 
region has a larger energy and axial angular momentum 
than a particle on a circular orbit in the region r > nsco > 
the region bound by the OSCO, if present, is likely to be 
vacant of particles on circular orbits (see discussion in 
Gair et al. 2008). Generally, however, this region is very 
close to the black hole, where the quasi-Kerr metric is no 
longer perturbative (see Paper I). 

The shaded region in Figure 2] marks the part of the 
parameter space where all circular equatorial orbits are 
radially stable. Albeit interesting in its own merit, we do 
not further consider this part of the parameter space, be- 
cause higher order terms in the radial epicyclic frequency 
become important. 

For particles in an accretion disk, circular motion at 
radii smaller than the ISCO is unstable, and the Kep- 
lerian frequency reaches its maximal value at the ISCO. 
The maximum Keplerian frequency is the largest oscilla- 
tion frequency in the disk and, therefore, an upper bound 
on the dynamical frequencies. In Figure |51 we plot the 
maximal Keplerian frequency ri0. max as a function of 
the spin a and the parameter e. Since the Keplerian 
frequency depends only weakly on the parameter e (see 
Figure [T|), quadrupole deviations affect the maximal Ke- 
plerian frequency predominantly due to the shift of the 
ISCO. The shaded region, again, marks the excluded part 
of the parameter space where higher order terms in the 
radial epicyclic frequency become important. 

3. QUASI-PERIODIC OSCILLATIONS 

In this section, we use our expressions of the Keplerian 
and epicyclic frequencies that we derived in Section 2 
to demonstrate how QPOs observed in black-hole X-ray 
binaries can be used to test the no-hair theorem in two 
different scenarios. In both interpretations of detected 
QPO signals, the dynamical frequencies are important 
but play different roles. We argue that either one or 
two of the three parameters mass, spin, and quadrupole 
moment can be determined from a pair of QPOs and that 
the no-hair theorem can be tested in conjunction with an 
independent measurement of the remaining parameters. 

3.1 THE DISKOSEISMOLOGY MODEL 

In general relativity, modes can be trapped in the ac- 
cretion disk of a black hole giving rise to modulations of 
the measured flux. Detailed expressions of these modes 
have been derived in a hydrodynamic model in full gen- 
eral relativity (the diskoseismology model; Perez et al. 
1997; Silbergleit et al. 2001). In this interpretation, the 
observed QPO pair can be identified as the lowest or- 
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der gravity (g-modes) and corrugation modes (c- modes) . 
The frequencies of these modes depend almost exclu- 
sively on the mass and spin in the case of a Kerr black 
hole and, consequently, both mass and spin can be deter- 
mined from a pair of QPOs (e.g., Wagoner et al. 2001). 

The radial epicyclic frequency K r has a maximum near 
the inner edge of the accretion disk (Okazaki et al. 1987). 
In Kerr geometry, the lowest order g-modes occur at fre- 
quencies that are very close to that maximum (Perez et 
al. 1997). Corrugation modes, on the other hand, can 
exist only in disks that corotate with the Kerr black hole 
with a spin a in the range < ciq < a < 0.95, where 
a ~ 10~ 5 - 10~ 3 (Silbergleit et al. 2001). The fun- 
damental c-mode coincides with the Lense-Thirring fre- 
quency of the black hole evaluated at a radius r c which 
is typically very close to the ISCO. For a black hole with 
a mass of 10M Q and a disk luminosity L ~ 10% of the 
Eddington luminosity, Silbergleit et al. (2001) estimate 
the radius r c by the formula 

r c = Hsco + K Q a- Kl {l-a) K \ (30) 

where K = 0.093M, K% = 0.79, and K 2 = 0.20. 

In a quasi-Kerr spacetime, the frequencies of the g- and 
c-modes depend not only on the mass and spin of the 
black hole, but also on the quadrupole moment through 
the parameter e. Therefore, the no-hair theorem can be 
tested if either the mass or spin is measured indepen- 
dently, as we demonstrate in the following. 

For our analysis, we make the simplifying assumptions 
that the frequency of the fundamental g-mode is equal 
to the maximum of K r /27r and that the fundamental c- 
mode is equal to the Lense-Thirring frequency of a quasi- 
Kerr black hole evaluated at the ISCO. For a Kerr black 
hole, the first assumption is an excellent approximation 
(Perez et al. 1997), while the second assumption provides 
an estimate that slightly deviates from the true value 
flur(r c )/2TT. In this working definition, we also assume 
that quasi-Kerr c-modes exist as long as Our (nsco) > 0. 
For a rigorous treatment, a detailed study of perturba- 
tions on accretion disks in quasi-Kerr spacetime is nec- 
essary, which is beyond the scope of this paper. 

In Figure [71 we plot as a function of the black hole spin 
(top panel) the fundamental frequencies of the g- and 
c-modes of a Kerr black hole using the above working 
definition (solid lines) as well as the fundamental c-mode 
(dashed line) evaluated at the radius r c given by equation 
(|30[) . The curves for the c-mode deviate primarily at 
lower spins. In Figure[7Jwe also plot (bottom panel) this 
deviation between the expressions for the fundamental 
c-mode as a function of the spin a. In the range of spins 
that is relevant for this paper, i.e., < a/M < 0.4, the 
deviation is ~ 10 — 30% except for very low values of the 
spin a < 0.1M. 

In Figure |S1 we plot (top panel) the fundamental g- 
mode and (bottom panel) the fundamental c-mode as a 
function of spin for several values of the parameter e in 
units of a 10M Q -black hole. In both cases, the frequency 
increases with increasing values of the spin, but decreases 
with increasing values of the parameter e. 

In Figure |H we plot frequency contours of the lowest 
order g-mode as a function of spin and the parameter e. 
The frequency of the fundamental g-mode increases for 
increasing values of the spin and decreasing values of the 
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Fig. 7. — Top: Fundamental frequency of the g-mode and c- 
modes for a Kerr black hole as a function of spin. We evaluate the 
fundamental c-mode at the radius of the ISCO (solid line) and at 
the radius r c (dashed line). Bottom: Percent deviation between 
both expressions of the fundamental c-mode as a function of spin. 



parameter e. The shaded region marks the excluded part 
of the parameter space. 

As can be seen from Figure O contours of constant 
frequency of the fundamental g-mode are fairly linear 
in the spin a and the parameter e. In order to further 
quantify this, we set r — > tq + Sa + rje, where ro = SM is 
the location of the maximum radial epicyclic frequency 
for a Schwarzschild black hole and S and r\ are constants. 
Then, we expand the derivative of the radial epicyclic 
frequency given by expression (|26[) with respect to radius 
around its maximum r = ro for a = e = and solve the 
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Fig. 8. — Frequencies of (top) the fundamental g-mode and (bot- 
tom) the fundamental c-mode scaled appropriately for a IOMq- 
mass black hole as a function of spin for several values of the pa- 
rameter e. Both the fundamental g-mode and c-modc frequencies 
increase with increasing values of the spin but decrease with in- 
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From these expressions, we derive a linear approximation 



Fig. 9. — Frequency contours of the fundamental g-mode as a 
function of the spin a and the parameter e scaled for a lOMQ-black 
hole. The fundamental g-modc frequency increases with increasing 
values of the spin and decreasing values of the parameter e. The 
shaded region is excluded from the parameter space. 
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Fig. 10. — Frequency contours of the fundamental c-mode as a 
function of the spin a and the parameter e scaled for a lOMQ-black 
hole. The fundamental c-mode frequency increases with increasing 
values of the spin and decreasing values of the parameter e. The 
shaded region is excluded from the parameter space. 

of the fundamental g-modc frequency in terms of the spin 
a and the parameter e given by the equation 
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where a* = a/AI . This approximation coincides with the 
maximum of the full expression of the radial epicyclic 
frequency given by equation (f!M| to within 10% across 
the full range of the values of the spin and the parameter 
e considered in Figure [9] 

In Figure 1 1 01 we plot frequency contours of the funda- 
mental c-modes as a function of spin and the parameter 
e. The c-mode frequency rises with increasing values of 
the spin and decreasing values of the parameter e. The 
shaded regions mark the excluded parts of the parameter 
space. The top left corner is excluded according to our 
requirement that quasi-Kerr c-modes only exist in coro- 
tating disks. Since the Lense-Thirring frequency given 
by equation (|28|) is the difference between the epicyclic 
frequencies in the <$> and 9 coordinates, its main depen- 
dence is on the parameter e. 
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Fig. 11. — Lines of constant g- and c-mode frequency (at 117.5Hz 
and 11.5Hz, respectively) as a function of mass and spin for several 
values of the parameter e. The intersection points of each corre- 
sponding pair of lines marks the particular combination of mass 
and spin for a given value of the parameter e. 

In Figures [TT] and [T^] we illustrate how measurements 
of g- and c-modes of a particular black hole can be used 
to test the no-hair theorem. For this purpose, we choose 
a black hole with a mass of M = 7M© and a spin of 
a = 0.2M. According to equations (J26]) and ([28]). the 
corresponding frequencies of the <?-mode and the c-mode, 
respectively, are 117.5Hz and 11.5Hz. In the following, 
we assume that these frequencies have been measured 
and identified as the respective lowest order g- and c- 
modes. 

In Figure 111! we plot contours of constant g- and c- 
modes for frequencies of 117.5Hz and 11.5Hz, respec- 
tively, as a function of the black hole mass M and spin 
a for several values of the parameter e. The intersection 
point of each pair of lines marks the particular combi- 
nation of the parameters M, a, and e that is consistent 
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Fig. 12. — Dependence of the black hole mass on the param eter 
e determined by the intersection points as illustrated in Figure [TT] 
Increasing values of the parameter e correspond to a larger mass. 
A separate measurement of the mass breaks the degeneracy of the 
mass and the parameter e allowing for a test of the no-hair theorem. 

with this hypothetical measurement. Increasing values 
of the parameter e shift the intersection point to larger 
values of the mass and spin. In the case of a Kerr black 
hole, this reproduces the assumed values for the mass 
and spin. 

Since we have three free parameters and two measure- 
ments, we can express the mass in terms of the param- 
eter e. In Figure [12J we plot the mass as a function of 
the parameter e for the intersection points given in Fig- 
ure 1111 Increasing values of the parameter e correspond 
to a larger black hole mass. In combination with an in- 
dependent mass measurement, the degeneracy between 
the mass and the parameter e is broken, and all three 
parameters mass, spin, and quadrupole moment of the 
black hole can be determined. This, in turn, allows us to 
test the no- hair theorem via relation (Q~]). 

3.2 THE KINEMATIC RESONANCE MODEL 

QPOs can also be modeled as resonances among 
the Keplerian and epicyclic frequencies (Kluzniak & 
Abramowicz 2001; Abramowicz et al. 2003). Several 
pairs of QPOs have been observed in galactic black holes 
at a ratio of w 3/2 (see Remillard & McClintock 2006). 
The assumption of either a 1:2 or a 1:3 resonance be- 
tween the Keplerian frequency and the radial epicyclic 
frequency for the observed QPO pair in the binary GRO 
J1655-40 has lead to a determination of the spin of this 
black hole (Abramowicz & Kluzniak 2001). 

We demonstrate the potential of this approach to test 
the no-hair theorem by considering an example similar to 
GRO J1655-40. We follow the analysis in Abramowicz & 
Kluzniak (2001), but we consider only the range of spins 
—0.4 < a/M < 0.4 in order to be consistent with our 
framework. For a 1:2 resonance between the Keplerian 
and radial epicyclic frequencies, we set Cl^ = 300 Hz and 
k x = 150 Hz. The resonance occurs at a radius r ros which 
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Fig. 13. — Spin versus mass of a black hole for a QPO pair at 
frequencies 300 Hz and 450 Hz in 1:2 resonance for several values 
of the parameter e. The mass increases for increasing values of the 
spin and decreasing values of the parameter e. 
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Fig. 14. — Spin vers us t he parameter e for the intersection points 
of the lines in Figure [T3l for different values of the mass M. The 
degeneracy between spin and the parameter e can be broken by an 
independent spin measurement. 

is a function of mass and spin for a Kerr black hole and 
of mass, spin, and the parameter e for a quasi-Kerr black 
hole. 

In Figure I13[ we plot the mass as a function of spin 
for and k x in a 1:2 resonance for several values of the 
parameter e. The mass increases for increasing values of 
the spin and decreasing values of the parameter e. If the 
mass of the central object is known from an independent 
measurement, then only certain combinations of the spin 



and the parameter e are consistent with such a measure- 
ment. For this illustration, we assume hypothetical mass 
measurements of M = 5.0, 5.5, 6.0Af Q , respectively. 

In Figure I14[ we plot the parameter e as a function of 
spin for the intersection points of the lines in Figure [13] 
with these hypothetical mass measurements. Higher val- 
ues of the spin correspond to larger values of the param- 
eter e for a given black- hole mass. The degeneracy be- 
tween spin and quadrupole moment can be broken by an 
independent spin measurement from, e.g., rclativistically 
broadened iron lines (e.g., Brenneman & Reynolds 2009), 
which in turn would uniquely determine the parameter e. 
In this QPO model, therefore, the no-hair theorem can 
be tested in conjunction with additional measurements 
of mass and spin. 

4. CONCLUSIONS 

In this series of papers, we are investigating a frame- 
work to test the no-hair theorem with observations of 
black holes in the electromagnetic spectrum (Paper I) 
based on a quasi-Kerr metric (Glampedakis & Babak 
2006). Since, according to the no-hair theorem, mass and 
spin are the only free parameters of a black-hole space- 
time, measurements of the mass, spin, and quadrupole 
moment can be used to test this fundamental property 
of black holes in general relativity. 

In this paper, we derived expressions of the Keplerian 
and epicyclic frequencies in the quasi-Kerr spacetime and 
explored their properties. We showed that, for moderate 
values of the spin, the Keplerian frequency depends al- 
most exclusively on the mass and spin of a given black 
hole, while the radial and vertical epicyclic frequencies 
depend significantly on the value of the parameter e near 
the inner edge of the accretion disk in addition to their 
dependence on the mass and spin. We determined the 
location of the ISCO for the ranges of the parameters 
-0.4 < a/M < 0.4 and -0.5 < e < 0.5 and showed 
that the Lense-Thirring frequency for a particle at the 
ISCO changes by a factor of ~ 5 for a black hole with 
spin a = 0.2M as the parameter e is varied from —0.2 
to +0.2. We discussed the emergence of two additional 
stable orbits inside the ISCO and possible consequences 
for astrophysical black holes. 

We demonstrated how this approach can be applied to 
QPOs observed in galactic stellar-mass black holes and 
AGN focusing on two different models. If a pair of QPOs 
is observed and identified as the fundamental g- and c- 
modes, respectively, the mass, spin, and quadrupole mo- 
ment of that black hole can be measured if the mass is 
known from dynamical observations. If the QPO pair is 
viewed as a nonlinear resonance between the Keplerian 
and epicyclic frequencies, the parameter e can be mea- 
sured together with independent measurements of the 
mass and spin. 

In addition to galactic black holes, Sgr A* is another 
prime target for tests of the no-hair theorem. High- 
resolution observations of stars in close orbit around 
Sgr A* for over a decade have lead to a precise mass 
and distance measurement of Sgr A* (Ghez et al. 2008; 
Gillessen et al. 2009). Recent VLBI observations re- 
solved Sgr A* on horizon scales (Doeleman et al. 2008) 
and pointed the way towards the first image of Sgr A* 
in the near future (Fish & Doeleman 2009). In Paper II, 
we showed that a ring of light will surround the image of 
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a black hole and that its shape depends directly on the 
mass, spin, quadrupole moment, and inclination of the 
black hole. In particular, we showed that the diameter 
depends predominantly on the mass, while the displace- 
ment and the asymmetry of this ring are proportional to 
the quantities asini and esin 3 / 2 i, respectively, where i 
is the inclination angle of the angular momentum of the 
black hole with respect to a distant observer. We noted 
that one additional observable is required to break the 
degeneracy of the displacement and the asymmetry with 
the inclination. 

Here we argue that a full test of the no-hair theorem 
for Sgr A* may be accomplished with a combination of 
VLBI imaging of the ring of light (or, more generally, 
the shadow; see Paper II) and the observation of quasi- 
periodic variability in the emission from Sgr A* using 



VLBI techniques (Doeleman et al. 2009). Such variabil- 
ity is thought to arise from density fluctuations orbit- 
ing around the center of mass. The orbital frequency of 
these inhomogeneities is the Keplerian frequency fi^,. In 
this paper, we showed that, for moderate values of the 
spin, the dependence of fl^ on the parameter e is negligi- 
ble. Consequently, if high-frequency VLBI observations 
of Sgr A* (Doeleman et al. 2009) can measure the Ke- 
plerian frequency of an orbiting hot spot and if Sgr A* 
is not spinning rapidly, they will measure the spin of 
Sgr A* irrespectively of the particular value of the pa- 
rameter e. Therefore, measurements of the displacement 
and the asymmetry of the ring of light (or the shadow) 
of Sgr A* will uniquely determine the inclination and the 
parameter e allowing us to test the no-hair theorem. 

We thank J. Schnittman and R. Wagoner for useful 
comments. This work was supported by the NSF CA- 
REER award NSF 0746549. 
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